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Abstract 

The amplitude for jp — > (pp is calculated in a pQCD kx - factorization approach. The total 
cross section for this process is compared with HERA data. Total cross section, as a function of 
photon-proton energy and photon virtuality, is calculated. We also discuss the ratio of clI^t and 
the dependence on the mass of the strange quark. The amplitude for ^p — > (pp is used to predict 
the cross section for exclusive photoproduction of <p meson in proton-proton collisions. Absorption 
effects are included. The results for RHIC, Tevatron and LHC energies are presented. 
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I. INTRODUCTION 



It was shown recently how to calculate the cross section for exclusive production of heavy 
quarkonia in the 7p — >■ Vp [l|, 0] as well as in pp — > pVp 0, 0] collisions within the k±- 
factorization approach. The same formalism was also used recently to calculate exclusive 
production of the Z° boson jBf. In the latter case the corresponding cross section is however 
very small and cannot be measured at present luminosities. 

In order to describe the photoproduction process in a strict pQCD framework, at least 
a hard scale must be present, which could be for example either a quark mass for heavy 
quarkonia, or a large photon virtuality for deep inelastic diffractive processes. The k±- 
factorization formalism described in [l|, [2| allows, however, a smooth continuation from the 
hard to soft regime, and therefore can be also used to model the photoproduction of light 
vector mesons. This derives from the fact that the fc^-factorization is a momentum space 
version of the color-dipole approach, where the soft region is modelled by the behaviour of 
the dipole cross section at large dipole sizes 0. 

The ZEUS collaboration at HERA has measured exclusive production of <fi mesons in the 
7P — > <pp reaction [tJ and in 7*p — > <pp reaction j8[. Here we wish to test the results of the 
approach against the HERA data as well as to make predictions for the pp — > p<pp reaction 
for the Tevatron and the pp — > p4>p reaction for RHIC and LHC. 



II. PHOTOPRODUCTION PROCESS jp 



W 




JF(ir, k) = dG(x, K)/dlogK 
FIG. 1: A sketch of the amplitude for exclusive photoproduction — > <f>p. 

The amplitude for the reaction is shown schematically in FigJTJ This amplitude has been 
derived for the first time in Ref.[l[,j2[. The imaginary part of the amplitude can be writen 
as: 

ZmM = w ^^ff™ [ d -^a s ^)H^^Kk 2 ) [ ^V(A V ,A 7 ) , (2.1) 
Ait 1 J k J z(l — z) 



2 



where /(Ay, A 7 ) have the form: 
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In the expressions above ip v (z, k) is the meson light-cone wave function and 3>i, $2 are given 
by (ref. Q): 
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r + k r — k r + A/2 r — A/2 

1 = ~ (r + k) 2 + e 2 ~ (r - k) 2 + e 2 + (r + A/2) 2 + e 2 + (r - A/2) 2 + e 2 ' ^ 2 ' 5) 

where: r = k + (z — \) A and e 2 = m 2 + z(l — z)Q 2 . 

The quantity J r (xi,X2, Ki, K2) is the off-diagonal unintegrated gluon distribution, which 
for small A 2 , within the diffraction cone can be approximated as: 

, J3 A 2 \ 

J r (x l ,x 2 ,k 1 ,k 2 ) = F(x eff ,K 2 ) exp ^ — - — J, (2.6) 

where a; e // = c s ( — ^ — J , c s — 0.41 and the forward unintegrated gluon distribution is taken 

from Ref. [9| where it was found in the analysis of the deep-inelastic scattering data. 

In the forward scattering limit, i.e. for A = azimuthal integrations can be performed 
analytically (see j4J), and one obtains the following representation for the imaginary part of 
the amplitude for the transverse photons: 
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QmM T (W,A 2 = 0,Q 2 ) = W 2 ^p^2 [ — ^- / Trdk'M^k 2 ) (2.7) 

47r Jo z \t ~ z ) Jo 

—a s (q 2 )T( Xeff , n 2 ) U (z, k 2 ) W (k 2 , k 2 ) + Ax(z, k 2 ) W^k 2 , k 2 ) 



where 
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^o(^ 2 , k 2 ) = t^—- 2 - -====== , (2.10) 

k + e 2 ^{k 2 - e 2 - k 2 ) 2 + Ae 2 k 2 

^(k\^ = 1-^(1+ (2.11) 



2fc 2 V v /(P-e 2 -K 2 ) 2 + 4e 2 A; 2 / ' 
The imaginary part of the amplitude for longitudinal photons is given as: 



where 
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In our calculation we choose the scale of the QCD running coupling constant at: 
q 2 = max{/€ 2 , k 2 + m 2 }. 

The full amplitude, at finite A 2 , for transverse and longitudinal polarization can be 
written as: 

M T (W, A 2 , Q 2 ) = (i + p) %mM T (W, A 2 = 0, Q 2 ) exp ^ B ^ A ^ ? (2 .15) 

M L (W, A 2 , Q 2 ) = {i + p) SmM L (W, A 2 = 0, Q 2 ) exp pggl^! ) ; ( 2 .16) 

where p is a ratio of real to imaginary part of the amplitude 

^teM (7 , d\og ksmM/W 2 ) 

P= =tan £ — ) = tan ( — Ajp ) . (2.17) 

H SmM V2 aiogiy 2 / V2 w ) y ' 

Above B(W) is a slope parameter which in general depends on the photon-proton center- 
of-mass energy and is parametrized in the present analysis as: 

B(W) = B + 2a' eff log {^), (2.18) 

with:S = 11 GeV~ 2 , a' eff = 0.25 GeV" 2 , W = 95 GeV (0). 

Assuming exponential dependence of the amplitude on A 2 (see Eq. ( 12.161) ) the total cross 
section for transverse and longitudinal polarization can be writen as: 
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1 + p 2 
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M L (W,A 2 = 0,Q< 



cx T (7p <pp) = 3m — nK ' tj79 '-^ . (2.20) 
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The full cross section is a sum of these two components: 

crtot{lP ->■ #) = o"t(7P -> #) + e c"l(7P -> #)• ( 2 -21) 

In HERA kinematics, e ~ 1. 

Let us collect now the main formulas involving the qq light-cone wave-function of the 
vector meson The electronic decay width depends on the model of the meson wave function 
as: 



T(V ^e + e~) = —^ -gl. (2.22) 
3M, 3 



In our calculation we use leading-order approximation, i.e. we neglect a possible NLO K- 
factor. The parameter gy can be expressed in terms of the <fi - meson wave function as 

a a 

a AT r rl 3 k 

9v = IT J Wf {M + mc,) Mk2) ■ (2 ' 23) 

Following [l[ 0] we use two types of the wave function, the Gaussian one, representing a 
standard harmonic-oscillator type quark model: 

^ ls (k 2 ) = C 1 e W (-^f) , (2.24) 



and a Coulomb wave function, representative of models in which the wave function has a 
long high-momentum tail: 



^ is (k 2 ) = -yL - ■ 2>2 . 2 . (2.25) 
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The parameters of the wave function are obtained by fitting to the decay width into e + e~ 
and imposing the normalization condition: 

1 = ^ r k 2 dk AMi, 2 ls {k 2 ) . (2.26) 
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III. EXCLUSIVE PHOTOPRODUCTION OF IN pp AND pp COLLISIONS 

The bare amplitude (when absorption effects are ignored) can be written schematically 
[3] as: 

M^( Pl ,p 2 ) = e 1 -^F KXl (p 1 ,t 1 )M rh2 ^vh 2 (s2,t 2 ,Q 2 1 ) (3.1) 
2 p 2 

+e2—- 1 J r x'x 2 (p 2 , h)M Th ^vhx{siM, Ql) ■ 

Z2 C2 
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FIG. 3: A sketch of the exclusive *yp —> (ftp amplitude with absorptive corrections. 



Because of the presence of the proton from factors only small Q\ and Q\ enter the amplitude 
for the hadronic process. This means that in practise one can put Q\ — Q\ = 

The full amplitude for the pp — > p<fip, reaction which includes elastic rescatterings reads: 

M( Pl ,p 2 ) = J ^S el (k)M^( Pl -k,p 2 + k) (3.2) 
= M^( Pl ,p 2 )-5M( Pl ,p 2 ), 

where 

S el (k) = (27r)V 2 )(fc) - ~T(fc) , T(k) = aZ(s) exp ( - l -B el k 2 ) . (3.3) 

In practical evaluations we take B e i = 14 GeV" 2 , of^ = 52 mb for the RHIC energy W = 
200 GeV, B el = 17 GeV" 2 , of£ = 76 mb [11] for the Tevatron energy W = 1.96 TeV and 
B el = 21 GeV" 2 , aZ = 100 mb for the LHC energy W = 14 TeV. 

The absorptive correction to the amplitude can be written as the following convolution 
of the bare amplitude (M*^) and the amplitude for elastic pp or pp scattering as: 

5M( Pl ,p 2 ) = J ^-T(k)M^( Pl -k,p 2 + k). (3.4) 

The differential cross section is expressed in terms of the amplitude M as 

da = ^-^| M | 2 dydtidt 2 dip . (3-5) 

This formula is used below to calculate several differential distributions in pp — > p<ftp and 
pp — > pcpp reactions. 
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IV. RESULTS 




FIG. 4: Total cross section for the photoproduction *yp — > (ftp process as a function of the photon- 
proton center-of-mass energy. In this calculation the Gaussian wave function is used. The curves 
are described in the text. 



Let us start review of our results with the jp — > (ftp reaction. In FigJH we show the 
total cross section for 'jp — > <pp as a function of photon-proton center of-mass energy W for 
Q 2 = GeV 2 . Our results are compared with the corresponding HERA data jzj. In order 
to describe experimental data we treat the mass of the strange quark as a free parameter, 
although it was fixed at m s = 0.37 GeV in Ref. Q. We show results for three different values 
of the strange quark mass. The red solid (upper) line is for ms = 0.37 GeV, blue (lower) 
line for ms = 0.50 GeV and the black line (which goes through the data points) for ms = 
0.45 GeV. We can see that the result for ms = 0.45 GeV gives the best description of the 
experimental data, so below we shall present results only for this value of the strange quark 
mass. Our results for Q 2 -dependence look consistent with the ones previously obtained by 
Ivanov in his phD thesis. [2[ 



J? 10 3 
c 




W = 75 GeV 1 


2 10 2 
-& 


- 




t 10 






Cl 

£ 1 






b -1 

10 








10 1 

Q 2 + M 2 (GeV 2 ) 




10 10 

Q 2 + M 2 (GeV 2 ) 



FIG. 5: Total cross section as a function of photon virtuality for two different wave functions: 
Gauss (left) and Coulomb (right) panel respectively. Here m s = 0.45 GeV. The solid line is for 
transverse and longitudinal cross section and the dashed line is only for transverse cross section. 



In Figj5] we show the total cross section as a function of photon virtuality for photon- 
proton center-of-mass energy W = 75 GeV. The data points are taken from Ref. [§[ . The solid 
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line is for the sum of longitudinal and transverse cross sections (<jl + <jt) and the dashed 
line is for transverse cross section {(Tt) alone. It is important that both these components are 
included in the calculation. In the left panel we show results for the Gaussian wave function 
and in the right panel for the Coulomb wave function. We can see that the Gaussian wave 
function much better describes the experimental data than the Coulomb wave function. 
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FIG. 6: The ratio cri/ar as a function of photon virtuality. The red solid line is for the Gaussian 
wave function and the black solid line is for the Coulomb wave function. Here m s = 0.45 GeV. 

In Figj6]we show the ratio of the longitudinal cross section (er^) to the transverse cross 
section {(Tt) as a function of photon virtuality, again for W = 75 GeV. The red solid line 
is for the Gaussian and the black line is for the Coulomb wave function. Our results are 
compared with the HERA data from Ref. 




FIG. 7: Rapidity spectrum of exclusive <f> production at RHIC (left), Tevatron (middle) and LHC 
(right) energies. The thin black lines are without, the thick red lines with absorption included. In 
this calculation Gaussian wave function was used and m s = 0.45 GeV. 

Let us come now to the presentation of our results for exclusive production of meson in 
hadronic reactions. In FigJT] we show distributions in rapidity for the pp — > p<pp (Tevatron) 
and pp — > pcpp (RHIC, LHC) reactions without (black thin solid) and with (grey thick 
solid) absorption effects. The absorption effects are rather small and will be not included in 
calculations of other observables. 

In FigJH] we show separate contributions of photon-pomeron and pomeron-photon fusion 
mechanisms in rapidity distribution of <fi in pp (middle panel) and pp (left and right panels) 
collisions. Here we have shown only results without absorptive corrections. The green line is 
for photon-pomeron, the blue line is for pomeron-photon fusion mechanism and the red most 
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FIG. 8: The photon-pomeron and pomeron-photon contributions to the rapidity distribution for 
the RHIC (left), Tevatron (middle) and LHC (right) energies. Here Gaussian wave function was 
used, m s = 0.45 GeV and absorptive corrections are not included. 



upper line is a sum of the two components. In rapidity distribution the two contributions 
add incoherently [3j. 
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FIG. 9: The transverse momentum distribution for the RHIC (left), Tevatron (middle) and LHC 
(right) energies. 

In Figj9]we show distribution in transverse momentum of the exclusively produced meson 
<p. In the left panel we present results for RHIC energy, in the middle panel for Tevatron 
energy and in the right panel for the LHC energy. 
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FIG. 10: Distribution in azimuthal angle for the RHIC (left), Tevatron (middle) and LHC (right) 
energies. 

Finally in FigflU] we show the distribution in relative azimuthal angle between outgoing 
protons. Quite different distributions are obtained for the RHIC, LHC (pp collisions) and 
Tevatron (pp collisions) . This effect is of the interference nature and was already discussed 
for the exclusive production of J/*f> jjjj]. 
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V. CONCLUSIONS 



In the present analysis we have performed the calculation of the cross section for the 
7P —7- <pp process as well as for exclusive production of <fi meson in the pp — > p<pp reaction. 

The cross section for the first process depends strongly on the model of the <ft meson 
wave function as well as on the strange quark mass. Both Gaussian and Coulomb wave 
functions were used with parameters adjusted to reproduce the electronic decay width of (f). 
The strange quark mass was adjusted to reproduce the jp —> <fip cross section measured by 
the ZEUS collaboration at HERA in the Gaussian model of the wave function. 

So-fixed parameters were used to make predictions for the pp — > pp(p and pp —> p<f>p 
reactions. We have presented distributions in the meson rapidity, in the meson transverse 
momentum and in the relative azimuthal angle between outgoing protons. The azimuthal 
angle dependence found here is of purely interference nature. Quite different distributions 
have been predicted for the Tevatron and LHC due to the different sign of the interference 
term (different electric charge of proton and antiproton). 

The cross sections found in the present analysis are fairly large. In practice one could 
measure the <fi meson either in the lepton-antilepton or in kaon-antikaon decay channels. The 
later method could be used by the ALICE collaboration at the LHC where even small trans- 
verse momenta of charged kaons can be registered A feasibility of the measurements 
requires, however, detailed Monte Carlo studies. 
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